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Abstract 

A geometric inequality in General Relativity relates quantities that 
have both a physical interpretation and a geometrical definition. It is 
well known that the parameters that characterize the Kerr-Ncwman 
black hole satisfy several important geometric inequalities. Remark- 
ably enough, some of these inequalities also hold for dynamical black 
holes. This kind of inequalities play an important role in the char- 
acterization of the gravitational collapse, they are closed related with 
the cosmic censorship conjecture. Axially symmetric black holes are 
the natural candidates to study these inequalities because the quasi- 
local angular momentum is well defined for them. We review recent 
results in this subject and we also describe the main ideas behind the 
proofs. Finally, a list of relevant open problem is presented. 



1 Introduction 



Geometric inequalities have an ancient history in Mathematics. A classical 
example is the isoperimetric inequality for closed plane curves given by 

> 47r^, (1) 



where A is the area enclosed by a curve C of length L, and where equality 
holds if and only if C is a circle (for a review on this subject see |9l]). General 
Relativity is a geometric theory, hence it is not surprising that geometric 
inequalities appear naturally in it. As we will see, many of these inequalities 
are similar in spirit as the isoperimetric inequality ([T]). However, General 
Relativity as a physical theory provides an important extra ingredient. It is 
often the case that the quantities involved have a clear physical interpretation 
and the expected behavior of the gravitational and matter fields often suggest 
geometric inequalities which can be highly non-trivial from the mathematical 
point of view. The interplay between geometry and physics gives to geometric 
inequalities in General Relativity their distinguished character. 

A prominent example is the positive mass theorem. The physics suggests 
that the mass of the spacetime (which is represented by a pure geometrical 
quantity [H] [IS] [30] ) should be positive and equal to zero if and only if the 
spacetime is fiat. From the geometrical mass definition, without the physical 
picture, it would be very hard to conjecture this inequality. In fact the proof 
turn out to be very subtle [T02] [TT8] . 

A key assumption in the positive mass theorem is that the matter fields 
should satisfy an energy condition. This condition is expected to hold for 
all physically realistic matter. It is remarkable that such a simple condition 
encompass a huge class of physical models and that it translates into a pure 
geometrical condition. This kind of general properties which do not depend 
very much on the details of the model are not easy to find for astrophysical 
objects (like stars or galaxies) which usually have a very complicated struc- 
ture. And hence it is difficult to obtain simple geometric inequalities among 
the parameters that characterize them. 

In contrast, black holes represent a unique class of very simple macro- 
scopic objects that play, in some sense, the role of 'elementary particles' in 
the theory. The black hole uniqueness theorem ensures that stationary black 
holes in electro- vacuum are characterized by three parameters, which can 
be taken to be the area A of the black hole, the angular momentum J and 
the charge q. The mass m is calculated in terms of these parameters by 
an explicit formula (cf. equation ([s])). It is well known that these parame- 
ters satisfy certain geometrical inequalities which restrict the range of them. 
These inequalities are direct consequences of the explicit formula ([s]). Among 
them, we note first the following 



which will lead to the Penrose inequality for dynamical black holes. Also 
we have the following two inequalities which will play a central role in this 




(2) 
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article 

^2 > q' + V'f + '^J' ^ ^ > An^/q^ + AJ^. (3) 

The equality in ^ is achieved for the Schwarzschild black hole. The equality 
in both inequalities ^ are achieved for extreme black holes. 

However black holes are not stationary in general. Astrophysical phenom- 
ena like the formation of a black hole by gravitational collapse or a binary 
black hole collision are highly dynamical. For such systems, the black hole 
can not be characterized by few parameters as in the stationary case. In 
fact, even stationary but non- vacuum black holes have a complicated struc- 
ture (for example black holes surrounded by a rotating ring of matter, see 
the numerical studies in [7]). Remarkably, inequalities (|2|-(|3]) extend (under 
appropriate assumptions) to the fully dynamical regime. Moreover, inequal- 
ities ([2])-(|3| are deeply connected with properties of the global evolution of 
Einstein equations, in particular with the cosmic censorship conjecture. The 
main subject of this review is to present a series of recent results which are 
mainly concerned with the dynamical versions of inequalities (|3]). 

To extend the validity of inequalities (|2])-([3]) to non-stationary black holes 
the first difficulty is how to define the physical parameters involved, most 
notably the angular momentum J of a dynamical black hole. To define 
quasi-local quantities is in general a difficult problem (see the review |1U8] ). 
However, for axially symmetric black holes, the angular momentum (via Ko- 
mar's formula) is well defined and it is conserved in vacuum. Essentially for 
this reason inequalities ^ have been mostly studied for axially symmetric 
black holes. An exception are inequalities which involves only the electric 
charge since the charge is well defined as a quasi-local quantity without any 
symmetry assumption. 

The plan of the article is the following. In section[2]we describe the heuris- 
tic physical arguments that support these inequalities and connect them 
with global properties of a gravitational collapse. In section [3] we present 
an overview of the main results concerning these inequalities that have been 
recently obtained. We also describe the main ideas behind the proofs. Two 
important geometrical quantities involved in ([3]), mass and angular momen- 
tum, have distinguished properties in axial symmetry (in contrast to the 
electric charge). These properties play a fundamental role. We describe in 
some detail the angular momentum in section |4] and the mass in section |5| 
Finally in section [6] we present the relevant open problem in this area. 
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2 The physical picture 



The most important example of a geometric inequality for dynamical black 
holes is the Penrose inequality. In a seminal article Penrose [95] proposed a 
physical argument that connects global properties of the gravitational col- 
lapse with geometric inequalities on the initial conditions. For a recent re- 
view about this inequality see [HE] and references therein. Since it will play 
an important role in what follows, let us review Penrose argument. 

We will assume that the following statements hold in a gravitational col- 
lapse: 

(i) Gravitational collapse results in a black hole (weak cosmic censorship). 

(ii) The spacetime settles down to a stationary final state. We will further 
assume that at some finite time all the matter have fallen into the black 
hole and hence the exterior region is electro-vacuum. 

Conjectures (i) and (ii) constitute the standard picture of the gravitational 
collapse. Relevant examples where this picture is confirmed (and where the 
role of angular momentum is analyzed) are the collapse of neutron stars 
studied numerically in [16] [58] . 

Before going into the Penrose argument, let us analyze the final stationary 
state postulated in (ii). The black hole uniqueness theorem implies that the 
final state is given by the Kerr-Newman black hole (we emphasize however 
that many important aspects of the black holes uniqueness still remain open, 
see [21] for a recent review on this problem). Let us denote by mo, Aq, 
Jo and go the mass, area, angular momentum and charge of the remainder 
Kerr-Newman black hole. In order to describe a black hole, the parameters 
of the Kerr-Newman family of solutions of Einstein equations should satisfy 
the remarkably inequality 

d>0, (4) 

where we have defined 

d = ml-ql- (5) 

'"'0 

This inequality is equivalent to 

From Newtonian considerations (take g = for simplicity), we can interpret 
this inequality as follows (see [109j ): in a collapse the gravitational attraction 
(~ ttIo/t^) at the horizon (r ~ mo) dominates over the centrifugal repulsive 
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forces ( ~ jQ/mor^). It is important to recall that the Kerr- Newman solution 
is well defined for any choice of the parameters, but it only represents a black 
hole if the inequality ^ is satisfied. 

The Kerr-Newman black hole is called extreme if the equality in ^ is 
satisfied, namely 

The area of the black hole horizon is given by the important formula 

Ao = 47r (2ml - ?o + 2moV^) . (8) 

Note that this expression has meaning only if the inequality ^ holds. 

Penrose argument runs as follows. Let us consider a gravitational collapse. 
Take a Cauchy surface S in the spacetime such that the collapse has already 
occurred. This is shown in figure [TJ Let S denotes the intersection of the 
event horizon with the Cauchy surface 5* and let A be its area. Let (m, q, J) 
be the total mass, charge and angular momentum at spacelike infinity. These 
quantities can be computed from the initial surface S. By the black hole area 
theorem we have that the area of the black hole increase with time and hence 

^0 > A. (9) 
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Since gravitational waves carry positive energy, the total mass of the space- 
time should be bigger than the final mass of the black hole 



m > mo- (10) 

The difference m — mo is the total amount of gravitational radiation emitted 
by the system. 

The area Aq of the remainder black hole is given by equation ([s]) in terms 
of the final parameters ( Jq, Qq, wiq). It is a monotonically increasing function 
of mo (for fixed go and Jo); namely the derivative dAQ/dmo is positive (we 
explicitly calculate this derivative bellow). Then, using this monotonicity 



and inequalities l\9h and (10) we obtain 



A<Ao<A-K\2m^-ql + 2m{m^-ql-^] 1, (11 



2 \ 1/2' 



m? 



where the important point is that in the right hand side appears the total 
mass m (instead of mo), which can be calculated on the Cauchy surface S. 
The parameters go and Jo are not known a priori but since they appear with 
negative sign we have 

2 \ 1/2' 



A < Ao < A-K \ 2m^ - ql + 2m (^m^ - ql - \^ j < lOvrm^. (12) 

Remains still an important point: how to estimate the area A of S in terms of 
geometrical quantities that can be locally computed on the initial conditions. 
Recall that in order to know the location of the event horizon the whole 
spacetime is needed. Assume that the surface S contains a future trapped 
two-surface Sq (the trapped condition is a local property). By a general result 
on black hole spacetimes we have that the surface Eo should be contained in 
S. But that does not necessarily means that the area of So is smaller than 
the area of S. Consider all surfaces S enclosing Sq. Denote by 74min(So) the 
infimum of the areas of all such surfaces. Then we clearly have that v4(Zl) > 
74min(So). The advantage of this construction is that y4min(So) is a quantity 
that can be computed from the Cauchy surface S. Using this inequality and 



inequalities (11) and (12) we finally obtain the Penrose inequality 



M^in(So) ^ (13) 

For further discussion we refer to [SB] and references therein. 

Penrose argument is remarkable because its end up in an inequality that 
can be written purely in terms of the initial conditions. On the other hand. 
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the proof of such inequahty gives indirect evidences of the vahdity of the 
conjectures (i) and (ii). 



Can we include the parameters q and J in the inequality (13) to get a 
stronger version of it? The problem is, of course, how to relate the final state 
parameters (go, Jo) with the initial state ones (g, J). 

If the matter fields are not charged, then the charge is conserved, namely 

q = qo- (14) 

And hence in that case we have the following version of the Penrose inequality 
with charge 

^mm(So) <A<An (^2m2 - q"^ + 2m (m^ - q^) ^^^j . (15) 

The case of angular momentum is more complicated. Angular momentum 
is in general non-conserved. There exists no simple relation between the total 
angular momentum J of the initial conditions and the angular momentum 
Jo of the final black hole. For example, a system can have J = initially, 
but collapse to a black hole with final angular momentum Jo ^ 0. We 
can imagine that on the initial conditions there are two parts with opposite 
angular momentum, one of them falls in to the black hole and the other scape 
to infinity. 

Axially symmetric vacuum spacetimes constitute a remarkable exception 
because the angular momentum is conserved. In that case we have 

J = Jo- (16) 



We discuss this conservation law in detail in section 5.1 The physical inter- 



pretation of (16) is that axially symmetric gravitational waves do not carry 
angular momentum. 

For non-vacuum axially symmetric spacetimes the angular momentum 
is no longer conserved. Matter can transfer angular momentum even in 
axial symmetry. This is also true for the electromagnetic field. However 
in the electro- vacuum case a remarkably effect occurs. First, the sum of 
the gravitational and the electromagnetic angular momentum is conserved. 
Second, at spacelike infinity only the gravitational angular momentum is 
non-zero. In section |4] we prove this two facts. Hence if J and Jq denotes 



now the total angular momentum, then the conservation (16) still holds for 
axially symmetric electro- vacuum spacetimes and we obtain the full Penrose 
inequality valid for axially symmetric electro- vacuum initial conditions 

^min(So) < A<47T \ 2m^ -q^ + 2m (m^ " " ^) 1 • (1^) 
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We emphasize that in this inequahty the total angular momentum J can 
be computed at any closed surface that surround the black hole using the 



formula (122). When this surface is at infinity the angular momentum is 
given by the gravitational angular momentum (i.e. the Komar integral of 
the axial Killing field). 



Inequality (17) implies the bound 



m > . (18) 



Of course, inequality (18) can be deduced directly with the same argument 
without using the area theorem. The first place where this conjecture was 
formulated is in |5l] (see also [73]). 



Inequality (18) can be viewed as a simplified version of the Penrose in- 



equality. The major difference is that the area of the horizon does not appears 



in (18). Only charges, which are essentially topological, appear in the right 



hand side of this inequality. 



Inequality (18) is a global inequality for two reasons. First, it involves 
the total mass m of the spacetime. Second it assumes global restrictions on 
the initial data: axial symmetry and electro-vacuum. We will discuss these 
assumptions in more detail in section |3j 

The area A, the angular momentum J in axial symmetry, and the charge q 



are quasi-local quantities (in particular, the right hand side of ( 18 ) is purely 
quasi-local). Namely they carry information on a bounded region of the 
spacetime. In contrast with a local quantity like a tensor field which de- 
pends on a point of the spacetime or a global quantities (like the total mass) 
which depends on the whole initial conditions. A natural question is whether 
dynamical black holes satisfy purely quasi-local inequalities. The relevance 
of this kind of inequalities is that they provide a much finer control on the 
dynamics of a black holes than the global versions. 

It is well known that the energy of the gravitational field can not be 
represented by a local quantity. The best one can hope is to obtain a quasi- 
local expression. These are the so called quasi-local mass definition (see the 
review article |108j and reference therein). Consider the formula ([s]) for the 
horizon area for the Kerr-Newman black holes. From this expression, we can 
write the mass in terms of the other parameters as follows 



In this equation we have dropped the subindice in the right hand side and 
also denoted the mass in the left hand side by rrihh to emphasize that this 
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expression can be in principle defined for any black hole (i.e. not necessar- 
ily stationary). With this interpretation, this expression is known as the 
Christodoulou I2HI mass of the black hole. 



For a dynamical black hole the expression (19) is in principle just a def- 



inition. Does the formula (19) represents the quasi-local mass of a non- 
stationary black hole? Let us analyze its physical behavior. We discuss first 
the relation of this quasi-local mass with the total mass of the spacetime. 
For only one black hole we expect the following inequality to be true 

m > nibh- (20) 



This inequality implies the Penrose inequality (17) but it is stronger (see the 
discussion in ^88j). However, it is important to emphasize that for the case 
of many black holes this inequality does not holds. In fact it is possible to 
find counter examples if we take the dbTQdb clS additive |117] or the quasi-local 
masses as additives [51]. This is expected, since the interaction energy of the 



black holes need to be taken into account (see the discussion in EH)- For 



only one black hole, the inequality (20) in axial symmetry is an open problem 



in this general form, we will further discuss it in section [6j 

We discuss now the purely quasi-local properties of ((T9|). The formula 



(19) trivially satisfies the inequality (18). This is, of course, just because the 



Kerr black hole satisfies this bound. Hence, if we accept (19) as the correct 



formula for the quasi-local mass of an axially symmetric black hole, then (19) 



provides the, rather trivial, quasi- local version of (18). 

Consider the evolution of rribh- By the area theorem, we know that the 
horizon area will increase. If we assume axial symmetry and electro- vacuum, 
then the total angular momentum (gravitational plus electromagnetic) will 
be conserved at the quasi-local level. On physical grounds, one would expect 
that in this situation the quasi-local mass of the black hole should increase 
with the area, since there is no mechanism at the classical level to extract 
mass from the black hole. In effect, the only way to extract mass from a black 
hole is by extracting angular momentum through a Penrose process. But 
angular momentum transfer is forbidden in electro-vacuum axial symmetry. 
Then, one would expect that both the area A and the quasi-local mass rribh 
should monotonically increase with time. 

Let us take a time derivative of rribh (denoted by a dot). To analyze this, 
it is illustrative to write down the complete differential, namely the first law 
of thermodynamics 

Srubh = ^SA + VIhSJ + ^nSq, (21) 
ovr 
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where 



4m 



bh 



An 
~A 



{q' + 4 J2 



AttJ 



H 



Am, 



H 



bh 



A 



(22) 

where rribh is given by (19) and d (defined in equation ([s])) is written in terms 
of A and J and q as 



d 



A 



(23) 



In equation (21) we have followed the standard notation for the formulation 



of the first law, we emphasize however that in our context this equation is a 



trivial consequence of ( 19 ) 



Under our assumptions, from the formula (19) we obtain 

rrihh = ^A, 



(24) 



were we have used that the angular momentum J and the charge q are 
conserved. Since, by the area theorem, we have 



A>0, 



(25) 



the time derivative of rrihh will be positive (and hence the mass nibh will 
increase with the area) if and only if /t > 0, that is 



A^T^/q^ + AJ^ < A. 



(26) 



Then, it is natural to conjecture that (26) should be satisfied for any black 



hole in an axially symmetry. If the horizon violate (26) then in the evolution 



the area will increase but the mass nibh will decrease. This will indicate 
that the quantity rubh has not the desired physical meaning. Also, a rigidity 



statement is expected. Namely, the equality in (26) is reached only by the 
extreme Kerr black hole given by the formula 



A = 471 + 4J2 



(27) 



The final picture is that the size of the black hole is bounded from bellow 
by the charge and angular momentum, and the minimal size is realized by 
the extreme Kerr- Newman black hole. This inequality provides a remarkable 
quasi-local measure of how far a dynamical black hole is from the extreme 
case, namely an 'extremality criteria' in the spirit of [21], although restricted 
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only to axial symmetry. In the article [13] it has been conjectured that, 
within axially symmetry, to prove the stability of a nearly extreme black 
hole is perhaps simpler than a Schwarzschild black hole. It is possible that 
this quasi-local extremality criteria will have relevant applications in this 



context. Note also that the inequality (26) allows to define, at least formally, 



the positive surface gravity density (or temperature) of a dynamical black 
hole by the formula (22) (see Refs. [15] [H] for a related discussion of the 



first law in dynamical horizons). 



If inequality (26) is true, then we have a non trivial monotonic quantity 
(in addition to the black hole area) rubh in electro-vacuum 



rhbh > 0. 



(28) 



It is important to emphasize that the physical arguments presented above 



in support of (26) are certainly weaker in comparison with the ones behind 
the Penrose inequalities (13), (15) and (17). A counter example of any of 



these inequality will prove that the standard picture of the gravitational col- 



lapse is wrong. On the other hand, a counter example of (26) will just prove 



that the quasi-local mass (19) is not appropriate to describe the evolution of a 



non-stationary black hole. One can imagine other expressions for quasi-local 
mass, may be more involved, in axial symmetry. On the contrary, reversing 



the argument, a proof of (26) will certainly suggest that the mass (19) has 



physical meaning for non-stationary black holes as a natural quasi-local mass 
(at least in axial symmetry). Also, the inequality (26) provide a non trivial 
control of the size of a black hole valid at any time. 

Finally, it is important to explore the physical scope of validity of these 
geometrical inequalities. Are they valid for other macroscopic objects? The 



Penrose inequality (17) is clearly not valid for an arbitrary region in the 



spacetime. Namely, consider an arbitrary 2-surface of area A, which is not 
necessarily a black hole boundary. We can make A arbitrary large keeping the 
mass (total or quasi-local) small (for example, take a region in Minkowski). 
For the inequalities (jTsj) and (26) the situation is less obvious. To have an 



intuitive idea of the order of magnitude involved it is important to include the 
relevant constants in these inequalities. Let G be the gravitational constant 
and c the speed of light. Then, these inequalities are written as follows 



and 



2 ^ 1 g2 + + 4 J2c2 

2 • 



G 



4vr^Vg4 + 4J2c2 < A. 



(29) 
(30) 
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For the reader's convenience we include the exphcit values of the constants 
(in centimeters, grams and seconds) 



G = 6.67 X 10"* cm-^g-'s', c = 3 x lO''^ cm s 



-1 



(31) 



The values of the fundamental physical constants used in this section were 
taken from [901. 



Let us analyze first the global inequality (29). It is useful to split it into 
the cases with zero charge and zero angular momentum respectively, namely 



Vg > 



m 

m 



and 



c ~ 



(32) 



(33) 



For the inequality (32) consider an electron and a proton, for these particles 
the quotient on the right hand side is given by 



rur. 



0.53 X 10^^ g-^^^cm^/'^s-^ 
2.87 X 10^^^ 



(34) 
(35) 



Since in these units we have vG ~ 10~ , we see these particles grossly violate 
the global inequality (32). And hence ordinary charged matter would violate 
it also (a similar discussion has been presented in [52] and [75j). 

For the angular momentum case (33) we can also consider an elementary 
particle. In that case the angular momentum of a particle with spin s (recall 
that s = 1/2 for the electron and the proton), is given by 



J = ^s(s + l)n, h = 1.05 X 10"^^ cm'^s-^g, 
where h is the Planck constant. For example, for the electron we have 

Ul 



Since 



G 

c 



1.12 X lO^^cm^s-^^"^ 



2.22 X 10~^^ cm^s'^g'\ 



(36) 



(37) 



(3^ 



inequality (33) is also violated by several order of magnitude for elementary 
particles. Instead of an elementary particle we can consider an ordinary 
rotating object. It is clear that there exists ordinary object for which J/m^ ~ 
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1 (say a rigid sphere of mass Ig, radius 1cm, and angular velocity Is ^) and 



hence inequality (33) is also violated for ordinary rotating bodies. 



We conclude that ordinary matter does not satisfies in general the global 



inequality (29 ). This inequality should be interpreted as a property of electro- 
vacuum gravitational fields on complete regular initial conditions where both 
the charge and the angular momentum are "produced by the topology" and 
not by matter sources (unless, of course, that they are inside a black hole 
horizon). By "produced by the topology" we mean the following. The an- 
gular momentum and the electric charge are defined as integral over closed 
two dimensional surfaces. In electro- vacuum these integrals are conserved 
(we discuss this in detail in section |4]) and hence they are zero if the topol- 
ogy of the initial conditions is trivial (i.e. M.^). In order to have non-trivial 
charges in electro- vacuum the initial conditions should have some "holes". 
This non-trivial topology signals the presence of a black hole. 



For the quasi-local inequality (30) it is also convenient to distinguish 



between the cases with zero charge and zero angular momentum respectively 

A > ^4V, (39) 



and 



A>8n%\J\. (40) 



C3' 



Since the charge is discrete, in unit of qe, it make sense to calculate the 
following characteristic radius 

a G^/"^ 

ro = ^ = 1.38 X 10"^^ cm. (41) 
c^ 

We see that tq is one order of magnitude less than the Planck length Ip given 

by 

lp={^^ =1.6x10-33 cm. (42) 

If we assume that the particle or the macroscopic object has spherical 
shape we can define the area radius r hy A = Attt"^. Then, the inequality 



(39) for a particle of charge Qe has the form 

r > ro. (43) 

The proton charge radius is Vp = 0.87 x 10"^^ cm according to [96] and 
r„ = 0.84 X 10"^^ cm according to the recent calculation presented in 



Hence, the proton satisfies inequality (43). This inequality is also consistent 
with the upper bound for the electron radius 10"^° cm measured in [52] . 
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For the case of angular momentum, using the relation (36) we can com- 
pute the following quotient for an elementary particle 



1/2 



^{s{s + l)f/%. 



(44) 



We see, that for a particle of spin s of order 1 the minimal radius is of the 
order of the Planck length Ip and hence it is also satisfied for elementary 
particles. 

More relevant is the case of an ordinary rotating body. The angular 
momentum J of a rigid body is given by 



J = lu, 



(45) 



where / is the moment of inertia and u the angular velocity. Consider an 
ellipsoid of revolution with semi-axes a and b, rotating along the b axis. The 
moment of inertia along the axis of rotation is given by 



(46) 



where m is the mass of the ellipsoid, which is assumed to have constant 
density. The area of the ellipsoid satisfies the following elementary inequality 



(47) 



The equality in (47) is achieved in the limit 6 — )■ 0, namely, the bound is 
sharp. 



Inserting equations (45) and (46) in to the inequality (40), and using the 



bound (47) and the fact that it is sharp, we obtain that the inequahty (40) 



is satisfied by the ellipsoid if and only if the following inequality holds 



5c3 
8G 



> mu. 



It is interesting to note that in the inequality (48) neither the area nor the 



radii of the body appear. The inequality relates only the mass and the angu- 
lar velocity of the body. Also, the body is not assumed be nearly spherical, 
the parameters for the ellipsoid are arbitrary. 

The value of the left hand side of this inequality is 



5c3 
8G 



2.53 X s'^g. 



(49) 
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For the sun we have the following values 



m,„„ = 1.989 X 10^3^, u = 2.90xl0-^rads'\ (50) 

and hence 

mu = 5.77 X 10"^^ s-^g. (51) 



We see that the inequality is satisfied for the sun. In order to violate (48) 
a body should be very massive and highly spinning, a natural candidate for 
that is a neutron star. For the fastest spinning neutron star found to date 
(see [To]) we have 

w ^ 4.5 X lOVarfs"^ (52) 

Assuming that the neutron stars has about three solar masses (which appears 
to be a reasonable upper bound for the mass, see [HS]) we obtain 



2.7 X lO^^s'i^. (53) 



The inequality (48) is still satisfied however the value (53) is remarkable close 



to the upper limit (49). 

The example of the ellipsoid shows the elementary relation between shape 
and angular momentum in classical mechanics which is valid even for non 
spherical bodies. There is no such relation between shape and the electric 
charge of a body. In fact, we will present some counter examples for charged 



highly prolated objects that violate the inequality (39). However, remarkably 



enough, for charged 'round surfaces' (we will define this concept later on), an 



inequality between area and charge can be proved (see theorem 3.5). On the 



other hand, the example of the ellipsoid suggests that the scope of validity 



of the inequality between area and angular momentum (40) (or the related 



inequality (48)) for axial symmetric bodies is much larger. 



3 Results and main ideas 



In this section we present the main results concerning inequalities ( 18 ) and 



(26) that haven been recently proved in the literature. We also discuss the 



general strategy of their proofs. 
3.1 Global inequality 

The proof of the inequality between total mass and charge (namely, setting 



J = in (18)), which is valid without any symmetry assumptions, has been 



known for some time. The first proof was provided in [62] and [60] using 
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Figure 2: Initial data with two asymptotically fiat ends. 



spinorial arguments similar to the Witten proof of the positive mass theorem 
|118j . See also [73]. A related inequality was proved in [HI] with similar 
techniques. In [T^] the proof was generalized to include low differentiable 
metrics. For this inequality an interesting rigidity result is expected: the 
equality holds if and only if the initial data are embed into the Majumdar- 
Papapetrou static spacetime (see [65] for a discussion of these spacetimes in 
relations with black holes). The rigidity statement has also been established, 
but with supplementary hypotheses, in ^62j and [2Z]. Very recently a new 
proof was provided in [HI] which removes all remaining hypothesis. Also in 
this article a new approach is presented. The strategy is to combine the Jang 
equation method with the spinorial proof of the positive mass theorem. 

The inclusion of angular momentum in axial symmetry (which is the main 
subject of this review) involves complete different techniques. In particular 
no spinorial proof of these inequalities are available so far (see however |119j 
where a related inequality is proved using spinors). The first proof of the 



global inequality (18) (with no electric charge) was provided in a series of 
articles [H], [IQ], [39] which end up in the global proof given in |44j . 

In [21] and [33] the result was generalized and the proof simplified. In 
|35j [38] the charge was included. As a sample of the most general result 
currently available we present the following theorem proved in [35] and [38] . 

Theorem 3.1. Consider an axially symmetric, electro-vacuum, asymptoti- 
cally flat and maximal initial data set with two asymptotics ends. Let m, J 
and q denote the total mass, angular momentum and charge respectively at 
one of the ends. Then, the following inequality holds 



^,^q' + ^q' + AJ^ 



m' > " ' ' . (54) 

For the precise definition, fall off conditions an assumptions on the electro- 



vacuum initial data we refer to [35] and [38]. For simplicity, in 5.1 we discuss 
in detail only the pure vacuum case. 
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Recall that for asymptotically fiat initial data the total mass m, the 
total charge q and the total angular momentum J (without any symmetry 
assumption) are well defined as integrals over two-spheres at infinity for a 



given asymptotic end. That is, all the quantities involved in (54) are well 
defined for generic asymptotically fiat data which are not necessarily axially 
symmetric. However the inequality does not hold without the symmetry 
assumption. General families of counter examples have been constructed in 
|74j for pure vacuum and complete manifolds. 

Under the hypothesis of this theorem (namely, electro-vacuum and axial 
symmetry) both the angular momentum and the electric charge are defined 
as conserved quasi-local integrals (we discuss this in detail in section |4]). In 
particular, if the topology of the manifold is trivial (i.e. M^), then these quan- 



tities are zero and hence theorem 3.1 reduces to the positive mass theorem. 



In order to have non-zero charge or angular momentum we need to allow 
non-trivial topologies, for example manifolds with two asymptotic ends as it 
is the case in theorem 3.1 (see figure [2]). An important initial data set that 



satisfies the hypothesis of the theorem is provided by an slice t = constant in 
the non-extreme Kerr-Newman black hole in the standard Boyer-Lindquist 
coordinates. 

This theorem has three main limitations: i) the initial data are assumed 
to be maximal, ii) there is no rigidity statement, iii) the data are assumed 
to have only two asymptotic ends. Let us discuss these points in more detail. 
The maximal condition plays a crucial role in the proof since it ensure a 
positive definite scalar curvature. A relevant open problem is how to remove 
this condition, we will discuss it in more detail in section |6l 



Extreme Kerr-Newman initial data, which reach the equality in (54), is 
not asymptotically fiat in both ends. The data have a cylindrical end and 
an asymptotically fiat end (see figure [3]). Hence these data is excluded in 



the hypothesis of theorem 3.1 In order to include the equality case we need 
to enlarge the class of data. An example is given by the following theorem 
proved in which includes the rigidity statement. 

Theorem 3.2. Consider a vacuum Brill initial data set such that they satisfy 
condition 2.5 in l4^ . Then inequality 



m>V|J|, (55) 



holds. Moreover, the equality in (55) holds if and only if the data are a slice 
of the extreme Kerr spacetime. 

The precise definition of the Brill class of data can be seen in [1^. The 
main advantage of these kind of data is that they encompass both class of 
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Figure 3: The cylindrical end on extreme Kerr black hole initial data. 



asymptotics: cylindrical and asymptotic flatness. We discuss this in section 



5.1 The condition 2.5 (see j44j for details) mentioned in this theorem im- 
plies that the initial data have non trivial angular momentum only at one 
end, however multiple extra ends with zero angular momentum are allowed. 
This condition involves also other restrictions which are technical. In a very 
recent work |105j these technical conditions have been removed and also an 
interesting new approach to the variational problem is presented. 



The inequality (54) (which in particular implies (55)) is expected to hold 
for manifolds with an arbitrary number of asymptotic ends, this general- 
ization is probably to most important open problem regarding this kind of 
inequalities (we discuss this in detail in the final section [6]). There exist, 
however, a very interesting partial result In order to describe it, we 

need to introduce the mass functional Ai, this functional is defined in sec- 
tion |5j It plays a major role in all the proofs, as it is explained in section 



3.3 This functional represents a lower bound for the mass. Moreover, the 
global minimum of this functional (under appropriate boundary conditions 
which preserve the angular momentum) is achieved by an harmonic map 



with prescribed singularities. As we will see in section |3.3 this is the main 



strategy in the proofs of all the previous theorem which are valid for two 
asymptotic ends. Remarkably enough in [33] the existence and uniqueness 
of this singular harmonic map has been proved also for manifolds with an 
arbitrary number of asymptotic ends. In this article the following theorem 
is proved. 

Theorem 3.3. Consider an axially symmetric, vacuum asymptotically flat 
and maximal initial data with N asymptotic ends. Denote by rrii, Ji (i = 
1, . . . N) the mass and angular momentum of the end i. Take an arbitrary 
end (say 1), then the mass at this end satisfies the inequality 

mi > MiJ2, . . . , Jn) (56) 

where Ai{J2, ■ ■ ■ , Jn) denotes the numerical value of the mass functional Ai 
evaluated at the corresponding harmonic map. 
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This theorem reduces the proof of the inequahty with multiples ends to 
compute the value of the mass functional on the corresponding harmonic 
map and verify the inequality 

M{J2,...,Jn) > V\J^\- (57) 
We further discuss this theorem in section 13.3.11 



Strong numerical evidences that inequality (57) holds for three asymp- 
totic ends has been provided in [48j. The numerical methods used in that 
article are related with the harmonic map structure of the equations. We 
will describe them in section 15. 2[ 

All the previous results assume complete manifolds without inner bound- 
aries. The inclusion of inner boundary is important to prove the Penrose 



inequality with angular momentum (17). Boundary conditions in relation 



with the mass functional Ai where studied in in order to prove a ver- 
sion of Penrose inequality in axial symmetry. In [36] inner boundaries were 
also included and an interesting new lower bound for the mass is obtained 
which depends only on the inner boundary. Finally, we mention that in [80] 



numerical evidences for the validity of the Penrose inequality (17) has been 
presented. 

3.2 Quasi-local inequalities 

Quasi-local inequalities between area and charge has been proved in [51] 
for stable minimal surfaces on time symmetric initial data. The following 
theorem proved in [17] generalize this result for generic dynamical black 
holes. 

Consider Einstein equations with cosmological constant A 

G,, = 87r(r^^/^ + T^,) - Ag^,, (58) 

where T^J^ is the electromagnetic energy-momentum tensor defined in terms 
of the electromagnetic field F^^i, by (85). The electric charge of an arbitrary 



closed, oriented, two-surface S embedded in the spacetime is defined by (86). 



Theorem 3.4. Given a closed marginally trapped surface S satisfying space- 
time stably outermost condition, in a spacetime which satisfies Einstein equa- 



tions (58) with non-negative cosmological constant A and such that the non- 
electromagnetic matter fields T^i, fulfill the dominant energy condition, the 
following inequality holds: 

A > ATcq"^, (59) 
where A and q are the area and the charge of S. 
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For the definition of marginally trapped surfaces (which is standard) and 
the stably condition see |2] [3] [79] |17] (see also [66] [98] )• This theorem is 
a completely quasi-local result that applies to general dynamical black holes 
without any symmetry assumption. It is also important to emphasize that the 
matter is not assumed to be uncharged, namely it is allowed that V ^F^^ ^ 
(which is equivalent to V^T^J^^ 7^ ). The only condition imposed in the 
non-electromagnetic matter field stress-energy tensor T^;^ is that it satisfies 
the dominant energy condition. In this theorem it is also possible to include 
the magnetic charge and Yang-Mills charges (see [37] and [77]). 



In |107] an interesting generalization of theorem 3.4 is presented in which 
the stability requirement is removed at the expense of introducing the prin- 
cipal eigenvalue of the stability operator and also the cosmological constant 
(with arbitrary sign) is added. 



At the end of section |2] we have observed that a variant of inequality (59) 
is expected to hold for ordinary macroscopic charged object (which are not 
necessarily black holes) at least if they are 'round enough'. In fact there exists 



an interesting and highly non-trivial counter example to (59) for macroscopic 
objects. This counter example was constructed by W. Bonnor in [22] (see 
also the discussion in [17]) and it can be summarized as follows: for any given 
positive number k, there exist static, isolated, non-singular bodies, satisfying 
the energy conditions, whose surface area A satisfies A < kq^ . The body is a 
highly prolated spheroid of electrically counterpoised dust. From the physical 
point of view we are saying that for an ordinary charged object (in contrast to 
a black hole) we need to control another parameter (the roundness) in order 
to obtain an inequality between area and charge. Remarkably enough it is 
possible to encode this intuition in the geometrical concept of isoperimetric 
surface: we say that a surface S is isoperimetric (or 'round') if among all 
surfaces that enclose the same volume as S does, S has the least area. Then, 
based on the results proved in [22] the following theorem was obtained in [17] 
for isoperimetric surfaces. 

Theorem 3.5. Consider an electro-vacuum, maximal initial data, with a 
non-negative cosmological constant. Assume that S is a stable isoperimetric 
sphere. Then 

A > —q\ (60) 
where q is the electric charge of S. 



Note that inequality ( 60 ) has a different coefficient as ( 59 ) . We recall that 
the notion of stable isoperimetric surface is very similar to the case of a stable 
minimal surface: the differential operator is identical, the only difference is 
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that the allowed test functions should integrate to zero on the surface, this is 
precisely the condition that the deformations preserve the volume (see [T7]). 



It is also possible to prove interesting variants of theorem 3.4 which are 
valid for generic surfaces (i.e. not necessarily trapped or minimal) but in 
order to obtain these results global assumptions on the initial data should 



be made (that is, in contrast to theorem 3.4, these are not a purely quasi- 
local results): the two surfaces are embedded on initial conditions that are 
complete, maximal and asymptotically flat. The non electromagnetic matter 
fields are assumed to be non charged and they should satisfy the dominant 
energy condition on the whole initial data (see Theorem 2.2 in |47j ) . 

As in the case of the global inequahty, quasi-local inequalities with angular 
momentum involve different techniques in comparison with the pure charged 
case. Their study star very recently. 



The quasi- local inequality with angular momentum and charge (26) was 
first conjectured to hold in stationary spacetimes in [10]. In that article the 
extreme limit of this inequality was analyzed and also numerical evidences for 
the validity in the stationary case was presented (using the numerical method 



and code developed in [7]). In a series of articles [HH] [SH] the inequality (26) 
was proved for stationary black holes. See also the review article [6j. 

It is important to emphasize that the stationary non-vacuum case is 
highly non-trivial. The physical situation is, for example, a black hole sur- 
rounded by a ring of matter (which do not touch the black hole). To illustrate 
the complexity of this case we mention that the Komar mass (which is only 
defined in the stationary case) can be negative for these black holes (for the 
Kerr-Newman black hole is always positive), see [H] [9]- It is interesting to 
mention that for this class of stationary spacetimes there exists a remarkable 
relation of the form (SvrJ)^ + (47rg^)^ = A'^A^ , where A'^ and A~ denote 
the areas of event and Cauchy horizon. This result have been proved in the 
following series of articles [H] [HZ] [1]. 

In the dynamical regime, the inequality was conjectured to hold in [45j 
based on the heuristic argument mentioned in section [2j In that article also 
the main relevant techniques for its proof were introduced, namely the mass 
functional on the surface and its connections with the area (we discuss this 
in section 5.3). A global proof (but with technical restrictions) was obtained 
in yj |j55j|. The first general and pure quasi- local result was proven in [5U] . 
where the relevant role of the stability condition for minimal surfaces was 
pointed out: 

Theorem 3.6. Consider an axisymmetric, vacuum and maximal initial data, 
with a non-negative cosmological constant. Assume that the initial data con- 
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tain an orientable closed stable minimal axially symmetric surface S. Then 



v4>87r|J|, (61) 
where A is the area and J the angular momentum of S. Moreover, if the 



equality in (61) holds then A = and the local geometry of the surface S is 



an extreme Kerr throat sphere. 

The extreme throat sphere geometry, with angular momentum J, was 
defined in |15] (see also [Jj and L50J). This surface capture the local geometry- 
near the horizon of an extreme Kerr black hole and it is defined as follows. 
The sphere is embedded in an initial data with intrinsic metric given by 

70 = 4j2e-"orf^2 + sin^ Odcj)'^ , (62) 



where ctq is given by (201). Moreover, the sphere must be totally geodesic. 



the twist potential evaluated at the surface must be given by ujq defined by 



(201) and the components of the second fundamental 

KijC = KijU^n' = Kijr]^r]' = 0, (63) 

must vanish at the surface. Here Kij denotes the second fundamental form 
of the initial data, ra* the unit normal vector to the surface and r/* the axial 
Killing field. Note that the functions cxo and Uq, which characterize the 
intrinsic and extrinsic geometry of the surface respectively, depend only on 
the angular momentum parameter J. The geometry of axially symmetric 



initial data set are described in detail in section SJ. In particular, the twist 
potential is determined by the second fundamental form Kij, using equation 



(175). The adapted coordinates system used in (62) is defined is section 4.1 
This theorem has two main restriction: the first one is the maximal con- 
dition. The second is vacuum. Remarkably enough, it is possible not only to 
avoid both restrictions but also to provide a pure spacetime proof (that is, 
no mention of a three-dimensional hypersurface) of this inequality, in which 
axisymmetry is only imposed on S. This generalization is proved in |79j : 

Theorem 3.7. Given an axisymmetric closed marginally trapped surface S 
satisfying the (axisymmetry- compatible) spacetime stably outermost condi- 
tion, in a spacetime with non-negative cosmological constant and fulfilling 
the dominant energy condition, it holds the inequality 

A>87r|J|, (64) 

where A and J are the area and (Komar) angular momentum of S. If equality 
holds, then S is a section of a non-expanding horizon with the geometry of 
extreme Kerr throat sphere. 
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The concept of non-expanding horizon is explained in ^9\, it essentially 
means that the shear vanished at the surface. 



It is important to note that the angular momentum that appears in (64) 
is the gravitational one (i.e. the Komar integral). The matter fields have 
also angular momentum and it can be transferred to the black hole, however 



the inequality (64) remains true even in that case. In fact this inequality 
is non-trivial for the Kerr-Newman black hole, we discuss this in detail in 
section HI 

In [99j it has been pointed out that there exists a connexion between the 



global inequalities described in section 3.1 and the quasi-local inequalities. 



This is obtained by linking the relevant mass functional Ai and (we 



discuss these mass functional in section 3.3) 



Inequality (64) has played an important role in the proofs of the non- 
existence of stationary two black holes configurations (see j93j (32j ) . 

Finally, we mention that there exists very interesting generalization of 



(64) to black holes in higher dimensions [72] . 



3.3 Main ideas 

In this section we present the main ideas behind the proofs of the global 



inequality (54) and the quasi- local inequahty (64). This section should be 
consider as a guide in which the technicalities are avoided. In section |4] and [5] 
we discuss in details the main relevant properties of the angular momentum 
and mass in axial symmetry which constitute the essential part of the proofs. 

3.3.1 Global inequalities 

The starting point in the proof of the inequality ([54j) (we consider the case 



q = for simplicity) is the formula for the total mass m given by Eq. (184). 
This formula represents the total mass as a positive definite integral over a 
maximal (here is where the maximal condition plays a crucial role) initial 
surface. This integral representation is a generalization of the Brill mass 
formula discovered in |26] (in section [s] we further discuss this formula and 
provide the relevant references). The formula holds in a particular coordinate 



system which is called isothermal (see lemma 5.1 ) 



The integrand in Eq. (184) has two kind of terms: dynamical and sta- 



tionary. The dynamical terms vanished for an stationary solution like Kerr. 



The stationary part lead to the relevant mass functional Ai defined by (185). 



This functional provides an obvious lower bound to the total mass, namely 

m>M{a,u). (65) 
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The mass functional Ai depends on two functions a and u. The function a 



is essentially the norm of the axial Killing vector (see equation (179)). The 
function uj is the twist potential of the Killing vector (see section |4|). These 
two functions can be freely prescribed on the initial data. This is an impor- 



tant and far from obvious property since the constraint equations ( 162 )-( 163 ) 
should be satisfied. This fact allows to formulate a variational principle for 
the functional Ai. The other important ingredient for this variational prin- 
ciple is the behavior of the angular momentum. The angular momentum is 
prescribed by the value of u at the axis (see section|4]). Then, if the variations 
of u vanishes at the axis the angular momentum will be preserved. With this 



two ingredients, it is hence possible to reduce the proof of the inequality (54) 
to a pure variational problem for the functional Ai{a,uj). 

The second step of the proof is to solve this variational problem. This 
is the most difficult part and also the most interesting since it reveals the 
geometric properties of the mass functional A^. Let us discuss the general 
strategy of the proof. 

Let ((Jo, Wo) be the corresponding functions obtained from the extreme 
Kerr initial data with angular momentum J, then we have that 

M{ao,uJo) = y^\. (66) 



The heuristic discussed in section |2j and Eq. (67) suggest that the following 
inequality holds 

Mia,uj)>y^l (67) 
for all {cr,u) such that u has the same value at the axis as the function uq. 



Moreover, the equality in (67) is reached if and only if a = do and oj = ojq. 
This is precisely the variational problem. Note the variational problem is 
formulated purely in terms of the functional A^(cr, w), without any reference 



to the constraint equations (162)-(163) 



The first evidence that this variational problem will have the expected 
solution is that the Euler-Lagrange equations of the functional Ai are equiv- 
alent to the stationary axially symmetric Einstein equations, in particular 
extreme Kerr satisfies these equation [H]. The second evidence (which is 
harder to prove) is that the second variation of Ai is positive definite evalu- 
ated at extreme Kerr [lO]. To prove this positivity property it is crucial to 
make contact with the harmonic map theory (in this case, trough the Carter 
identity). With these ingredients it is possibly to show that extreme Kerr 



is a local minimum of the mass and hence the inequality (67) is proved in 
an appropriate defined neighborhood of extreme Kerr. This local proof was 
done in BOl. 
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To have a global proof of (67) (i.e. without any smallness assumption) 
more subtle properties of the mass functional are required. A crucial step is 
to realize that the mass functional Ai is essentially the renormalized energy 
of an harmonic map into the hyperbolic plane (we discuss this in section 



5.2). This kind of harmonic maps have been extensive studied in the litera- 



ture. The problem here is that the map is singular at the axis and hence the 
standard techniques do not apply directly. To use the harmonic map theory 
we need to handle these singularities and this is the main technical difficulty. 
In [H] the proof of this variational result was done using estimates, inspired 
in the work of G. Weinstein [m] [m] [Il3] [HI] [Il6] [Il5] (see also [56])), 
which rely on particular properties of this functional (inversion symmetry). 
In subsequent works [31], [33], [35], [38] the proof was simplified and im- 
proved using general results on harmonic maps (more precisely the existence 
result [7T]). In these proofs the connection with the harmonic maps theory 
is more transparent and the problem of the singularities is clearly isolated. 

Finally let us mention the following important point regarding initial data 
with multiples ends. The multiples ends appears in the variational problem 
(67) as singular points of the functions (cr, u). Remarkably, even in that case 
it is possible to solve completely the variational problem. This is precisely 



the content of theorem 3.3 (proved in [33]). The only missing piece in order 
to prove the inequality (54) in that case is the following. Theorem 3.3 ensure 



the existence of a global minimum but the value of Ai at the global minimum 
is unknown. In the case of two asymptotic ends we known that extreme Kerr 



is the global minimum and hence we can explicitly compute the value (66) 



3.3.2 Quasi-local inequalities 



The global inequality (54) applies to a complete three dimensional manifold. 



In contrast the quasi-local inequality (64) apply to a closed two-surface. In 



principle it is a priori not clear at all that there is a relation between these two 
kind of inequalities. The two physical heuristic argument presented in section 
|2]in support for them are very different. In particular, it is far from obvious 
that the mass functional Ai can play a role for the quasi-local inequalities. 
Remarkably enough, it turns out that a suitable adapted mass functional 
Ai'^ over a two-surface play a very similar role as Ai. The motivation for the 
definition of Ai'^ given by (204) is discussed in section 5.3, The new mass 



functional Ai'^ and its connection with the area represent the kernel of the 
proof of (64). Let us discuss this. 



On the extreme Kerr initial data there exist an important canonical two- 
surface, namely the intersection of the Cauchy surface with the horizon. On 
a surface given hj t = constant in the Boyer-Lindquist coordinates this two- 
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surface is located at infinity on the cylindrical end (see figure |3]). Let call 
(ctojCUo) the value of the functions (ctqjWo) (defined in the previous section 



3.3.1) on this two-surface. The functions {ao,uJo) will play for JH'^ a similar 
role as the functions {ao,uo) for Ai. Namely, first they satisfy the Euler- 
Lagrange equations of . Second, the second variation of is positive 
definite evaluated at ((10,^0) (see [IS]). There is also a connection with the 



harmonic maps energy (we discuss this in section 5.3). Using similar kind of 



arguments as described in section 3.3.1 it is possible to prove the following 
inequality 

2|J| <e(-^'^('^''^)-«)/^ (68) 

for all functions {a, u) such that u has the same values at the poles of the 
two-surface as loq. With equality if and only if we have a = ao and u = uq. A 



local version of this inequality was first proved in [4:5] . The global version (68 ) 



was proved in [Ij. The rigidity statement was proved in [5U]. We emphasize 



that the inequality (68) (in complete analogy to the inequahty (67)) is a 
property of the functional M'^ , the geometry of the two-surface S does not 
intervene at all. 



The inequality (68) is interesting but, in the light of the discussion pre- 
sented in 3.3.1 it is somehow expected. What is crucial and completely 
unexpected is the relation between M.^ and the area of the surface S. This 
relation was founded locally in fl5] and globally in [Tj. By local in this con- 
text we mean that the relation holds for surfaces in an appropriated defined 
neighborhood of an extreme Kerr throat geometry. On the other hand, global 
means that the relation holds for general surfaces. In jHO] the important con- 
nection with the stability condition was proved for minimal surfaces and in 
[79j for marginally trapped surfaces. The final inequality essentially reads as 
follows. For a two-surface S such that: (i) it is minimal [50] (or marginally 



trapped [79]) and (ii) it is stable, then the following inequality holds 

A > 47re(^'-«)/^ (69) 

The minimal and marginally trapped conditions are requirements on the 
extrinsic curvature of S (namely, on the trace of the second fundamental 
form). The stability condition is a requirement on derivatives of the second 
fundamental form. In the spacetime version presented in [79]) the inequality 



(69) is a consequence of a flux inequality (see Lemma 1 in that reference) 



where the geometric and physical meaning of each term is apparent. 



26 



4 Angular momentum in axial symmetry 



Axial symmetry plays a major role in the inequalities that include angular 
momentum presented in the previous sections for two main reasons. For 
the global inequality ( 54 ) is the conservation of angular momentum implied 
by axial symmetry which is relevant. For the quasi-local inequality (61), is 
the very definition of quasi-local angular momentum (only possible in axial 
symmetry) which is important. These two properties are closed related for 
vacuum spacetimes, since the Komar integral provides both the conservation 
law and the definition of quasi-local angular momentum. For non-vacuum 
spacetimes in axial symmetry the Komar integral still provides a meaningful 
expression for the gravitational angular momentum but it is no longer con- 
served. Nevertheless, as we already mentioned in section |2| in the electro- 
vacuum case the sum of the gravitational and electromagnetic angular is 
conserved, and hence in that case is possible to prove the global inequality 
(54). In the general non- vacuum case in axial symmetry (that is, with general 
matter sources which are not electromagnetic) no simple and universal rela- 
tion between mass and angular momentum is expected. For example in 
neutron starts models in axial symmetry have been numerically constructed 
such that they violate inequality (54). Nevertheless, remarkably, the quasi- 
local inequality (61) for black holes is still valid in the general non- vacuum 
axially symmetric case. 

In this section we summarize relevant results concerning angular momen- 
tum in axial symmetry. Although these results are not new they are not 
easy to find in the literature, notably, the conservation of angular momen- 
tum in the electro-vacuum case and the relation between Komar integrals 
and potentials in the presence of matter fields. The conservation of angular 
momentum in the electro-vacuum, together with the relevant Komar inte- 
gral for the electromagnetic angular momentum, were discovered in [106j . in 
a much general setting. We present a simpler derivation of this result in the 
language of differential forms. 

Let us begin with some general remarks about conserved quantities in 
General Relativity (see also the review article [TS] for a related discussion and 
|lU8j for the general problem of how to define quasi-local angular momentum 
without symmetries). 

Let M be a four dimensional manifold with metric (7^1, (with signature 
( — h ++)) and Levi-Civita connection V^. On this curved background, let 
us consider an arbitrary energy-momentum tensor T^^ which satisfies the 
conservation equation 

V^T^"' = 0. (70) 
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It is well known that if the spacetime admit a Killing vector field r]^ 



V(^r/,) = 0, (71) 

then the vector 

= (72) 

is divergence free 

V^.K^' = 0. (73) 

This equation provides an integral conservation law via Gauss theorem. For 
some of the computations in this section it is convenient to use differential 
forms instead of tensors. We will denote them with boldface. Let K be the 



1-form defined by (72). Equation (73) is equivalent to 



d*K = 0, (74) 

where d is the exterior derivative and the dual of a p form is defined with 
respect to the volume element e^jyA7 of the metric Qf^i, by the standard formula 

Let Vt denotes a four-dimensional, orientable, region in M and let dVL be its 



three-dimensional boundary. Then using (74) and the Stokes theorem we 
obtain 



0= / d*K= / *K. (76) 
Jn Jan 

Note that in this equation the region Q is arbitrary and the boundary dQ 
can have many disconnected components. 

Consider a spacelike three-surface S. The conserved quantity correspond- 
ing to the Killing vector rj'^ is defined with respect to S by 

K{S) = I *K. (77) 
Js 



The interpretation of equation (76) in relation to the quantity K{S) is the 
following. Let f2 to be a timelike cylinder, such that its boundary dfl is 
formed by the bottom and the top spacelike surfaces 5*1 and 5*2 and the 
timelike piece C. Then we have (taking the corresponding orientation) 

0= f *K = K{Si)-K{S2)+ [*K. (78) 
Jan Jc 



The integral over the timelike surface C is the fiux of K{S). Equation (78) is 



interpreted as the conservation law of the quantity K{S). The region fl can 
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also be chosen to have a null boundary C, equation (78) remains identical 
and the interpretation is similar. 

If the Killing vector t]^ is also a symmetry of the tensor T^,^ (we have not 
assumed that so far), namely 

£rjTf.u = 0, (79) 

where £ denote Lie derivatives, then the following vector is also divergence 
free 

K, = M{T,X - ^Tt],), (80) 

where T is the trace of T^,y. In fact there is a whole family of divergence free 
tensors since Trj^ is divergence free. Hence the previous discussion applies 
to this vector as well. 

Note that the conserved quantities -^'(5') are naturally defined as inte- 
grals over spacehke three-surfaces. In flat spacetime it is possible to convert 
these integrals into a boundary integral over two-surfaces. This gives the 
quasi-local representation of conserved quantities (see the discussion in the 
introduction of |108j ). In a curved background this is in general not pos- 
sible. However, as we will see, this possible for the particular case of the 
electromagnetic field. 

Before analyzing the angular momentum it is important to study the elec- 
tric charge. The electric charge is of course relevant in our discussion since 
it appears in the inequalities discussed in the previous sections. But more 
important, even if we want to analyze these inequalities in pure vacuum, the 
electric charge represent the simpler 'conserved charge' on a curved space- 
time. Its definition and properties serve as model for all the other conserved 
quantities, like the angular momentum. 

The Maxwell equations on (M, g^^) are given by 

V^F^. = -47rj„ (81) 
V[^F,„] = 0. (82) 

In terms of forms, they are written as 

d*F = 47r*j, (83) 
dF = 0. (84) 

The energy-momentum tensor of the electromagnetic field is given by 

T^u = ^ (^F^xF/ - ^g,,Fx^F'^^ . (85) 
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Let 5 be a closed orientable two-surface embedded in M (in the following, 
all two-surfaces will be assumed to be closed and orientable). The electric 
charge g of 5 is defined by 



1 

An 



Let S" be a three-surface with boundary iS, then using Stokes theorem and 



Maxwell equation (83) we obtain 



q{S) 



(87) 



This equation is interpreted as follows. From equation (83) we deduce the 



conservation law for the current j analog to (74), namely 

= 0. (88) 
And hence taking the same region VL and using Stokes theorem we obtain the 



analog expression as ( 78 ) for the current j 



J+ / J- 

52 Jc 



19) 



Using (87) we finally obtain 



q{Si) - q{S2) 



J- 



(90) 



This is the conservation law for the electric charge. Note that in the left 



hand side of (90) we have integrals over two-surfaces, in contrast with (78) 



where integrals over three-surfaces appear. This is because we have an extra 



equation (i.e. Maxwell equation (83)) that allow us to write the integral (87) 
in the form (1861). 



(91) 



When j = the charge has the same value, namely 

q{Si) = q{S2), 



and we say that the charge is strictly conserved. 

We turn now to angular momentum for axially symmetric spacetimes. 
We begin with the definition of axial symmetry. 

Definition 4.1. The spacetime {M^g^y) is said to he axially symmetric if its 
group of isometries has a subgroup isomorphic to 5*0(2). 
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We will denote by rj'^ the Killing field generator of the axial symmetry. 
The orbits of rj'^ are either points or circles. The set of point orbits F is called 
the axis of symmetry. Assuming that F is a surface, it can be proved that rj'^ 
is spacelike in a neighborhood of F (see |H9]). We will further assume that 
the Killing vector is always spacelike outside F. Note that if this condition is 
not satisfied then the spacetime will have closed causal curved, in particular 
it can not be globally hyperbolic. 

The form rj^ will be denoted by r/, and the square of its norm by rj, namely 

V = V^V^ = \r]\^- (92) 

We have used the notation rj^ to denote the Killing vector field and rj to 
denote the square of its norm to be consistent with the literature. However, 
in this section, to avoid confusions between rj'^ and its square norm rj, we 
will denote the vector field 77^ by fj in equations involving differential forms 
in the index free notation. 

Consider now Einstein equations on an axially symmetric spacetime 

G^u = SttT^,. (93) 



Note that the Killing equation (71) implies that T^,^ satisfies (79). 

The Komar integral (it is also appropriate to call it the Komar charge) 
of the Killing field is defined over a two-dimensional surface S as follows 



(94) 



Hence, as we discussed above, via Stokes theorem we obtain 



J(S) = — I d*dr], (95) 



s 



167r 

where 5* is a three-dimensional surface with boundary S. It is a classical 



result [82] (see also PJDJ) that the integrand in (95) can be computed in 
terms of the Ricci tensor 

{e.a]/Sj'^^V^) = lR^.v''^^.a0y. (96) 
In terms of forms this equation is written as 

d*dr] = 2*K, (97) 
where we have defined the 1-form K by 

K = K^ = R^^r]". (98) 
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Using Einstein equations (93) the form K can be written in terms of the 
energy momentum tensor 



(99) 



Note that this expression is identical to (80). Then, repeating the same 



argument, we obtain the conservation law for angular momentum in axial 



symmetry which is the exact analog to the charge conservation (90) 



J(5i) - J{S2] 



Stt 



K. 



(100) 



The right hand side of this equation represent the change in the angular 
momentum of the gravitational field which is produced by the left hand side, 
namely the angular momentum of the matter fields. Note that the angular 
momentum of the matter fields are written as integrals over three-dimensional 
surfaces. In particular in vacuum we have the strict conservation of angular 
momentum 

J{Si) = J{S2). (101) 
We consider now the case where the energy-momentum tensor in Einstein 



equations (93 ) is given purely by the electromagnetic field (85 ). For simplicity 



we consider the case with no currents j = 0, which is the relevant one since 



only in that case we get conserved quantities. In that case we have that (85) 



satisfies the equation (70) and the source-free Maxwell equations are given 
by 

d*F = 0, (102) 
dF = 0. (103) 

We also assume that the Maxwell fields are axially symmetric, namely 

£r,F = 0. (104) 

Consider the 1-forms defined by 

cx = F{n), f3 = *F{rj), (105) 

where we have used the standard notation F{f]) = Fn^rj'^ to denote contrac- 
tions of forms with vector fields. Using the general expression for the action 
of the Lie derivative on forms 



£fjUJ = d[uj{r])] + {duj){v), 



(106) 
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Maxwell equations (102)-(103) and the condition (104) we obtain 



dct = 0, dl3 = 0. 
It follows that there exist locally functions x ^i-nd ip such that 

cx = dx, f3 = dip. 



(107) 



(108) 



The form K defined by (99) has the following expression for the electro- 
magnetic field 

1 



K = 2(F{a)--rt\F\' 



We have that (see |116] ) 

*(77A(F(a)) =aA/3. 



Using (108) we obtain 



d{a A /3) = 0, 
and hence there exist locally a 1-form 7 such that 

dj = cx A f3, 

where 7 is given by 



Note that 



7 = ^ (X# - i^dx) . 



From equation (110) we deduce 

*{r]AK) = 2ctA(3. 
We use the following identity valid for arbitrary 1-forms 
*{r]AK)A'n = r]*K - *r]{K{f])). 



Using (112) we finally obtain our main formula 

*K = 2d(7 A :^) + 27 A + *fi{K ■ 77), 



where we have defined 



It is important to note that 



rj 



dfi{r]) = 0. 



(109) 

(110) 
(111) 

(112) 

(113) 
(114) 

(115) 

(116) 

(117) 
(118) 
(119) 
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Using equation (117), we integrate *K over a three-surface 5* tangential to 
r/^, with boundary S. Using that r/^ is tangential to 5* it follows that the 



restriction of the 3-form *'r) to 5* is zero. For the second term in (117) we use 



equations (114) and (119) to obtain the same conclusion. Hence, we have 



2d{j A ff) 



2 ^V, 



(120) 



where in the last equality we have used Stokes theorem. 

We summarize the previous calculation in the following lemma, which is 
a re- writing of the result that have been obtained in |106] . 

Lemma 4.2. Consider an axially symmetric spacetime for which the Einstein- 



Maxwell equations (93), (85), (102), (105) are satisfied. Let S he an ori- 
entable three- surface, tangent to the axial Killing field rj^ , with boundary 
(possible disconnected) S. Then we have 



J{S) 



1 



7 Af), 



:i2ii 



where J{S) is the Komar integral given by (94), 7 is defined in terms of the 



electromagnetic field by (113) and r) is given by (118) 



We can define a 'total angular momentum' which is conserved in electro- 
vacuum, namely 

MS) = ^ [*dr]-4-fAfi. (122) 



We note that since the surface S is tangent to rj^ and all the fields are axially 
symmetric, then the surface integrals are in fact line integrals on the quotient 
manifold M\ 50(2). 

for rotating isolated horizons. This 



Formula (122) was studied in 



formula has been also recently studied at null infinity in connection with the 
center of mass and general definition of angular momentum for asymptoti- 
cally flat (at null infinity) spacetimes, see |83] . 



A very important example where lemma 4.2 applies is the Kerr-Newman 
black hole. Consider the Kerr-Newman black hole with parameters (m, a, q). 
The total angular momentum is given by Jt = am. This is equal to the 
Komar integral evaluated at infinity, since the electromagnetic field decay 
and does not contribute at infinity. However at the horizon the Komar an- 
gular momentum is not Jt- The angular momentum at the horizon has the 



decomposition (122). 
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For the Kerr-Newman black hole the Komar angular momentum at the 
horizon J is given by (see, for example, [9^ page 222). 



J = a 



2r+ \ 'Za^ \ ar_|_ 
where r+ is the horizon radius 

r_|_ = m + (m^ 
The area of the horizon is given by 

A = 47r(r^ + 



arctan 



2 2\l/2 

a — q ) ' . 



(123) 



(124) 



(125) 



tioned, it is well known that the area satisfy the inequality 



Equation (|125|) is of course identical to equation (|8|). As we already men- 

(126) 



A > 47rA/g4 + 4j2 



Which in particular implies 



A > SttIJtI. 



(127) 



But this inequality relates the total angular momentum J^- It is a priori not 
obvious if the following inequahty holds 



A > 87c\J\, 



:i28) 



where J is given by (123), namely the Komar integral at the horizon. Note 



that this is precisely the inequality proved in theorem 3.6 and that the Kerr- 



Newman black hole satisfies all the hypothesis of that theorem. 



Let us check explicitly that indeed (128) is satisfied. The expression 



(|123|) is remarkably complicated, to better analyze it let us rewrite it in the 

(129) 



following form. Using (125) we have 

J = 

where 



A 



in 



a 



-I- 



arctan 



(130) 



Instead of using a it is convenient to use x = a/r^ as free parameter. In 
terms of x the function e is written us 



X 1 + 



2r 



(131) 
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where 



X 



+ X arctanfx 



(132) 



We take as free parameters (g,r+,x). Note that — 1 < x < 1 and hence we 
have 



|e| < 



1 + 



2r 



(133) 



Fix (g,r+). It can be exphcitly check that function f{x) is non-positive and 
have a unique global minimum at x = where /(O) = —2/3. Hence it follows 
that 

|e| < 1. (134) 



4.1 Potentials 

The potentials for the axial Killing vector plays an important role in the 
mass functional described in section |5l 

It is instructive to analyze first the electric charge and its potential in 
axial symmetry. Assume first that the Maxwell equations are source free. 



Then we have found the potentials x and ip defined by equation (108). In 
particular, the potential ip determine the electric charge over an axially sym- 
metric two-surface S. In order to see that it is convenient to consider a tetrad 
(P, k^, 1]'^) and coordinate system (6*, 0) adapted to an axially symmetric 
two-surface defined as follows (see figure |4]). For simplicity we will assume 
that S has the topology of a two-sphere. Let us consider null vectors i'^ and 
fc^ spanning the normal plane to S and normalized as i'^k^^ = —1, leaving a 
(boost) rescaling freedom i'^ = fi'^, fc'^ = f~^k^. By assumption rj'^ is tan- 
gent to S, it has on the surface closed integral curves and vanishes exactly 
at two points which are the intersection of the axis F with S. We normalize 
vector T]'^ so that its integral curves have an affine length of 27r. Let us chose 
a coordinate on iS such that t]^ = d/dcj). The other vector of the tetrad 
which is tangent to S and orthogonal to rj^^ will be denoted by and assume 
that it has unit norm. We define the coordinate 9 such that is proportional 
to d/dO and such that = tt, are the poles of S. 

The induced metric and the volume element on S (written as spacetime 
projectors) are given by q^y = g^u + ^^lku + Lk^ and e^y = 2~'^ex^f,J^k^ 
respectively. The area measure on S is denoted by ds. Since the surface is 
axially symmetric we have Cj^q^y = 0. 

Using this tetrad, the charge over an axially symmetric surface S is writ- 
ten as follows 

g(5) = ^ / "^''^'^'f^^ds. (135) 
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Figure 4: Adapted tetrad for an axially symmetric two-sphere 



where is defined by (105). In the source free case we can use the potential 
ip defined by ( 108 ) to obtain 



'tt Js Jo 



'e^rf^. =^(7r) -V^(O), 



(136) 



That is, the charge is given by the difference of the value of the potential ip 
at the poles of the surface S. 

We have seen that the potential ip is only defined in the source-free case 
as scalar function in the spacetime. However, on the surface S it is always 
possible to define a potential ip by the equation 



(137) 



on the surface, since this equation involves only a derivative with respect 
to 9. The potential ip is only defined on the surface. By definition, in the 
source-free case where the potential ip also exists, the two functions are equal 
up to a constants on the surface. This constant is irrelevant since it does not 
affect the charge. 

Consider now the potential for the angular momentum. The twist vector 
of Tj^^ is defined by 

It is a well known result that the vacuum equations 



duj = 0, 



(138) 
imply that 

(139) 



where uj is the 1-form defined by (138). Hence there exist a potential u such 
that 

LJ = du. (140) 

The function u is the twist potential of the Killing vector r]^, it contains all 
the information of the angular momentum as we will see. In the non- vacuum 
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case, the twist potential is not defined. In the following we will not assume 
vacuum. 

In terms of the adapted tetrad we have that the Komar expression is 
given by 



J 



We the relation 



1 



1 



in 



Vrj'^l^K ds. 



A., 7 



-1„ 



to write the Komar integral in the following form 

1 



J 



IGtt 



It is clear that the vector defined by 

^7 = V''^'^,uxyi''k\ 



:i4ii 



(142) 



(143) 



(144) 



is orthogonal to l^^k^^ r]^ and it has unit norm and hence is the other member 
of the tetrad. Then we have 

J=^ [ r"'i'^,ds. (145) 
ISvr Js 

We emphasize that this expression is valid only for axially symmetric surfaces. 

If we assume vacuum, then we can use the twist potential defined by ( 140 ) 
to obtain 



J 



deoodO 



(a;(7r)-a;(0)). 



(146) 



Equations (145) and (146) are the analog to equations (135) and ( 136[ ) for 
the charge. 

We consider now the non-vacuum case. We define the following vector 
(a part of the extrinsic curvature of S with the role of a connection on its 
normal cotangent bundle, see e.g. the discussion in [63] [23] ) 



Since S is axially symmetric, the tetrad can be chosen such that 



= 0. 



In particular this implies that 



(147) 



(148) 



(149) 
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Using this equation we obtain that 



(150) 



From this equation we get the following equivalent expression for the Komar 
angular momentum (see, for example, 



J 



1 

8^ 



^fiW ds. 



(151) 



Remarkably, as we will see in the following, the vector fi^f^ determines even 
in the non-vacuum case a potential on the surface (we will denote it by u) 
which coincide with the twist potential u defined above in the vacuum case. 
We emphasize that the potential u will be only defined on the two-surface 
S, in contrast to the twist potential u which (in the vacuum case) is defined 



by equations (139) and (140) on a region of the spacetime. 



By construction, the vector fi^f'' is tangent to the surface S. Since we have 
assumed that S has the topology, there exists functions u and A such that 
the vector has the following decomposition on S in terms of a divergence-free 
and an exact form 



n^j^ = Ta + VaA, 



Ta 



(152) 



The functions u and A are fixed up to a constant. In equation (152) we have 



used the capital indices (which run from 1 to 2) to emphasizes that this is 
an intrinsic equation on the two-surface S. 

By assumption both and the intrinsic metric on S are axially sym- 
metric, then it follows that the functions u and A are also axially symmetric 
(i.e. they depend only on 6). In particular it follows that 



0, 



0, 



(153) 



where is the vector tangent to S previously defined. Since the norm rj is 



also an axially symmetric function, equation (153) implies that 

T^'DaV = 0. 



(154) 



function uj such that 



Equation (|154|) ensures the integrability conditions for the existence of a 

1 

2t]' 



T, 



A 



-eAB^^CO. 



(155) 



Note that equation (|155|) is valid only in axial symmetry. Collecting these 

(156) 



results we obtain the decomposition 



n 



27] 



eAsV^u + VaX. 
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In particular, using (153), we obtain 



1 

2r] 



(157) 



We write equation (157) with spacetime indices and we use the following 
representation for the tetrad vector 



= V 



to finally obtain 



Using equation (151) an integrating we finally get 

1 



J 



^(0)). 



(158) 



(159) 



(160) 



The relevance of this construction is that the function w, which is only defined 
at the surface iS, plays the role of a potential that can be defined in the non- 
vacuum case [79]. To see the relation between u and note that we have 



(161) 



This equation is valid always (i.e. non-vacuum) and it gives the function uj 
in terms of the twist vector w^. In the non-vacuum case the twist potential 
uj is not defined but the function u is always well defined. In fact, we can 



take (161) as definition of w, since the right hand side is only a derivative 



with respect to 9. On the other hand, in the vacuum case equation (161) 



implies that u and u differs by a constant which is irrelevant since it does 
not contribute to the angular momentum. 



5 Mass in axial symmetry 

The main goal of this section is to present the mass formula for axially 
symmetric data (184) and the mass functional for two-surfaces (204). We 



also discuss their geometrical properties in connection with harmonic maps. 



5.1 Axially symmetric initial data 



The global geometrical inequality (54) is studied on axially symmetric, asymp- 
totically fiat initial data set. In order to present the mass formula, we first 
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review the basic definitions and properties of tfiis kind of initial data. For 
simplicity we concentrate only in the pure vacuum case (see [35] [38] for the 
electro- vacuum case). 

An initial data set for the Einstein vacuum equations is given by a triplet 
{S, hij,Kij) where is a connected three-dimensional manifold, hij a Rieman- 
nian metric, and Kij a symmetric tensor field on S. The fields are assumed 
to satisfy the vacuum constraint equations 

DjK'^ - D'K = 0, (162) 
R - KijK'^ + = 0^ (163) 

where D and R are the Levi-Civita connection and the scalar curvature 
associated with h^, and K = Kijh^K In these equations the indices are 
moved with the metric hij and its inverse h^^ . 
The initial data is called maximal if 



K = 0. (164) 
This condition is crucial because it implies via the Hamiltonian constrain 



(163) that the scalar curvature R is non-negative. In fact on the right hand 



side of equation (163) it is possible to add a non-negative matter density and 



all the arguments behind the proof of theorem 3.1 will also apply since they 
rely on lower bounds for R. 

The manifold S is called Euclidean at infinity, if there exists a compact 
subset /C of 5* such that 5* \ /C is the disjoint union of a finite number of open 
sets Uk, and each Uk is diffeomorphic to the exterior of a ball in M^. Each 
open set Uk is called an end of S. Consider one end U and the canonical 
coordinates in which contains the exterior of the ball to which U is 
diffeomorphic to. Set r = (^(x*)^) . An initial data set is called asymptot- 
ically flat if S is Euclidean at infinity, the metric hij tends to the euclidean 
metric and K^j tends to zero as r — oo in an appropriate way. These fall off 
conditions (see \L8l [30] for the optimal fall off rates) imply the existence of 
the total mass m (or ADM mass [ill) defined at each end U by 

1 f 

m = —— lim / {djhij — dihjj) s^ds, (165) 

167r r->oo J 

where d denotes partial derivatives with respect to x\ Sr is the euclidean 
sphere r = constant in U, s* is its exterior unit normal and ds is the surface 
element with respect to the euclidean metric. 

The angular momentum is also defined as a surface integral at infinity 
(supplementary fall off conditions must be imposed, see, for example, the 
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review articles |108j . [78] and reference therein). Let /3* be an infinitesimal 
generator for rotations with respect to the fiat metric associated with the 
end U, then the angular momentum J in the direction of /3* is given by 

J = lim / (Ki,: - KhiAl3's^ds. (166) 

StT r-!>oo J 

The above discussion applies to general asymptotically flat initial data. 
We have seen that at any end the total mass m and the total angular mo- 
mentum J are well defined quantities. We consider now axially symmetric 
initial data. In analog way as the spacetime definition |4.1 , we say that the 



Riemannian manifold {S, hij) is axially symmetric if its group of isometries 
has a subgroup isomorphic to SO (2). We will denote by ?7* the Killing field 
generator of the axial symmetry and by F the axis. The initial data set is 
called axially symmetric if in addition 77* is also a symmetry of Kij, namely 

£nKij = 0. (167) 

On an axially symmetric spacetime there exist axially symmetric initial con- 
dition, conversely axially symmetric initial data evolve into an axially sym- 
metric spacetime. However, on an axially symmetric spacetime it is also 
possible to take initial conditions which are not axially symmetric in the 
sense defined above. These conditions will of course also evolve into an axi- 
ally symmetric spacetime, but the Killing vector is 'hidden' on them (these 
kind of initial data were studied in |2D])- We will not consider this kind of 
data here since on an axially symmetric spacetime it is always possible to 
chose initial conditions which are explicitly axially symmetric in the sense 
defined above. 

For axially symmetric data we have the Komar integral discussed in sec- 
tion |4j This integral can be calculated in terms on the initial data if we chose 
S G S. There exist a very simple relation expression for the Komar integral 
on an initial data, namely 

J = / Kij7]'s^ds. (168) 
0^ Js 



The equivalence between (168) and the original definition (94) can be seen 
as follows. Consider the tetrad adapted to S defined in section 4.1 Assume 
that S (Z S. By assumption, the axial Killing vector is tangent to the 
three-dimensional surface S. Denote by the timelike unit normal to the 
three-surface S", and let s'^ be the spacelike unit normal to the two-surface S 
which lies on S. These vectors can be written in term of ^^ and as follows 

n^=-L(£M + p)^ = _L(£M_fcM)^ (169) 
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Using these expression, we can write the integrand in (94) as 



J{S) = ^ / n\s^VxV,- 



(170) 



The second fundamental form can be written (as spacetime tensor) in terms 
of the unit normal as 



K 



Using that 'q^n^ = and that is tangent to S, from (171) we obtain 



:i7i] 



(172) 



From ( [172| ) and ( [170| ) we obtain (jTeSj). 

Comparing expression ( |168 ) with (166) we see that if we chose in (166) 
the vector /3* = r/* (that it, the generator of axial rotations) then these to 
expression are equivalent since rj^Si = for an sphere at infinity. 

It is possible to calculate the potential u for the spacetime Killing field 
defined in section |4] in terms of Kij as follows. Define the vector Ki by 



— ^iikS'^Tj , Si 



(173) 



where €ijk is the volume element with respect to the metric hij. Then, as a 

(174) 



consequence of equations (162) and (167) we have that 





Then the potential u is defined by 



Ki = --DiU. 



(175) 



It can be checked that this is the same potential as defined in the previous 
section (see [12]) evaluated on the initial surface. The importance of the 
expression ( |175 ) is that allow to calculate to in terms of the initial conditions. 

The functions u and rj have the information of the stationary part of the 
initial conditions. The dynamical part is characterized by the functions rj' 
and u' defined by 



V 



(176) 



The notation comes from the fact that they are related with the time deriva- 
tive of this functions, namely if n'^ is the timelike normal of the initial data, 
we have (see [42j) 

ry' = n'^V^T], u' = n^V^o;. (177) 
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The whole initial data can be constructed in terms of {ri,uj]ri' ,uj'). These 
functions are the initial data for the wave map equations that characterize 
the evolution in axial symmetry. This is clearly seen in the quotient repre- 
sentation (see P9]). For our present purpose the relevant part of the initial 
conditions is contained in {r],oj). 

So far we have discussed local implications of axial symmetry. Suppose 
that S is simply connected and asymptotically flat (with, possible, multiple 
ends). It can be proved (see [31]) that in such case the analysis essentially 
reduces to consider manifold of the form S = ^^\Ylk=o where ik are points 
in M'^. These points represent the extra asymptotic ends of 5*. Moreover, in 
[31] it has been proved that on 5* there exists the following global coordinate 
systems which will be essential in what follows. 

Lemma 5.1 (Isothermal coordinates). Consider an axially symmetric, asymp- 
totically flat (with possible multiple ends) Riemannian manifold {S, hij), where 
S is assumed to be simply connected. Then, there exist a global coordinate 
system system (p, z, ip) such the metric has the following form 



where the functions a,q,Ap,Az do not depend on if. In these coordinates, 
the axial Killing vector is given by rj = 8/ dip and the square of its norm is 
given by 



Using this coordinates system, the end points ik are located at the axis 
p = of M^. Define the intervals /fc, 1 < /c < — 1, to be the open sets 
in the axis between ik and i^-i, we also define Jq and Jat as 2 < iq and 
z > i^ respectively (see figure |5| . The manifold S is Euclidean at infinity 
with N + 1 ends. In effect, for each i^ take a small open ball of radius 
r(fc), centered at ik, where r(fc) is small enough such that does not contain 
any other ik' with k' ^ k. Let Br be a closed ball, with large radius R, such 
that Br contains all points ik in its interior. The compact set /C is given by 
/C = Br \ XlfcLi Bk and the open sets Uk are given by B^ \ ik, for 1 < k < N, 
and Uq is given by \ Br. Our choice of coordinate makes an artificial 
distinction between the end Uq (which represent r — )■ oo) and the other ones. 
This is convenient for our purpose because we want to work always at one 
fixed end. We emphasize that the isothermal coordinates (p, z, (p) (and hence 
the corresponding spherical radius r = y/p^ + z'^) are globally defined on S. 
In what follows we will always use these coordinates. 

The smoothness of the initial data on the axis implies that the potential u 
is constant on each interval Ik- These constants are directly related with the 



h = e("+2«)(rfp2 + dz^) + p^e^dif + Apdp + AJzf, 



(178) 




(179) 
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«0' 



Figure 5: Axially symmetric data with asymptotic ends. 



angular momentum of the end points ifc- In effect, the angular momentum 
of an end ifc is defined to be the Komar integral with respect to a surface 



that encloses only that point. Using the formula (146) we obtain 



J, ^ J(5,) = - ( 



-J- 



:i8o) 



Let (So a surface that enclose all the end points z^, then the total angular 
momentum of the end r — )■ oo is given by 



Jo = J{So) = I (a;|/„ - uj\i^) , 



(181) 



which is equivalent to 



N 



Jq — Jk 



(182) 



k=l 



The mass (165) is defined as a boundary integral at infinity. For axially 



symmetric initial data there is an equivalent representation of the mass as a 
positive definite integral over the three-dimensional slice. This formula was 
discovered by Brill [26] and allow him to prove the first version of the positive 
mass for pure vacuum axially symmetric gravitational waves. This formula 
was generalized to include metrics with non hypersurface orthogonal Killing 
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vectors (i.e. with non-zero Ap and Az in the representation (178)) in [6T] and 
for non-trivial topologies in jUj. In particular, this includes the topology 
of the Kerr initial data. See also |12] and [IH] for related discussion on this 
mass formula. 

Take isothermal coordinates and assume that the initial data are maximal 
condition ( 164[ ) holds). Then the total mass (165) is given by the 



I.e. 



following positive definite integral 



m 



3271 



\da\^ + 



:i83) 



0"^ + a\ and dv is the flat volume element in M^, namely 



where \da\'^ 

dv = pdzdpdif. Since the functions are axially symmetric the integration in 
is trivial, we keep however this notation because it will be important in 
section 15.21 



The remarkable fact of the integral (183) is that it applies also to the 



case of multiples ends, the integral is taking over all the extra ends since 
it cover the whole M^. 

The integrand in (183) can be further decomposed as follows (see [19] for 
details) 



327r 



{v 



/2 



u'^) +F^ + \da\ 



\du\ 

jj2 



dv. 



:i84) 



Where the term is essentially the square of the second fundamental form 
of the two-surface that is obtained by the quotient of the manifold S by the 
symmetry 5*0(2) (this two-surface is a half plane p > minus the end points 
ifc). The important point in this representation is that it clearly separates 
the dynamical and the stationary terms in the mass formula. In fact, this 
formula can be seen as the geometrical conserved energy of the wave map 
that appear in the evolution equations (see [39]). 



The first three terms in the integrand of (184) correspond to the dynam- 



ical part for the mass. They vanished for stationary solutions. The last two 
terms correspond to the stationary part. They give the total mass for Kerr. 
Explicitly the stationary part is given by 



327r 



|(9wp) dv, 



An hence we obtain the important bound 

m > M(a,u). 



:i85) 



;i86) 
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We have presented the mass formula but we have not discussed the asymp- 
totic conditions that ensures that this integral is well defined. We discuss 
them now trying to focus only in the essential aspects, for the technical de- 
tails we refer to and [31]. All the important features can be seen already 
in the term \da\'^ in (185). We will concentrate in this term only in what 
follows. 

It is instructive to see how the function a behaves for the Kerr black 
hole initial data with mass m. These initial data have two asymptotic ends, 
namely we have only one extra end point ii. In the limit r — > oo (that is, 
the end Uq) we have 

- + 0(r-2). (187) 



At the end we have 



a = -41og(r) + 0(1) 



for the non-extreme case, and 



a = -21og(r) + 0(l), 



:i89) 



for the extreme case. We see that the fall off behavior at the end ii is different 
for the non-extreme and extreme case, this reflect of course the difference 
between an asymptotically flat end and a cylindrical end (see flgures |2] and 
^respectively). 

To include this kind of fall off behavior let us consider the following con- 
ditions. In the limit r — )■ oo we impose 



a 



'1/2 



o(r-i/^), da = o(r-^/^). 



3/2> 



At the end point ik we impose 



da = o(r 



-3/2x 
(fc) J 



(190) 



(191) 



Despite the formal similarity of (191) and ( |190 ) they are of a very different 
nature. Conditions (190) are essentially the standard asymptotic flat con- 



ditions on the metric, they include in particular the behavior (187) of the 
Kerr initial data. On the other hand, conditions (191) on the extra ends 
ik are much more general than asymptotic flatness. Note that asymptotic 
flatness at the ends ik in these coordinates essentially implies a behavior like 
(190). But conditions (188) includes also (189) (in fact also any logarithmic 



behavior with any coefficient). 

The motivation for the fall off conditions (191) and (190) is that they 
ensures that the integral of |c?(Tp is finite both at infinity and at the end points 
ik- And the advantage of condition (191) is that it encompasses both kind 



of behavior: asymptotically fiat and cylindrical. This is the main motivation 



of the Brill data definition that appears in theorem 3.2 (see [44] for details). 
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5.2 Harmonic maps 



The crucial property of the mass functional defined in ( 185 ) is its relation to 
the energy of harmonic maps from to the hyperbolic plane H^: they differ 
by a boundary term. To see this relation consider first the mass functional 
A4 defined on a bounded region Q such that Q does not intersect the axis F 
(given by p = 0) 



327r 



{\da\' 



p-^e-^'^ldu]^) dv. 



(192) 



We consider this integral for general functions a and u which are not nec- 
essarily axially symmetric. The function log p is harmonic outside the axis, 
that is 

Alogp = 0, inM3\r, (193) 

where A is the flat Laplacian in M^. Using equation (193) we obtain the 
following identity 

d log p 



Mo = M' 



-(logp + a) ds, 



(194) 



where the derivative is taken in the normal direction to the boundary dQ 
and is given by 



M'o 



327r 



|9r/|2 + \du\ 



dv. 



(195) 



Recall that rj is defined by (179). 

The functional M'q defines an energy for maps {t], w) : — )■ where 
denotes the hyperbolic plane {{r],u) : r] > 0}, equipped with the negative 

constant curvature metric 



ds' 



drf + duj'^ 

^2 



The Euler-Lagrange equations for the energy Ai'^ are given by 

\du\^ 



A log f] = - 
Aw = 2 



^1 
dudr] 
V 



(196) 

(197) 
(198) 



The solutions of (197)-(198), i.e, the critical points of are called har- 
monic maps from M.^ — )• H^. Since A4n and A4q differ only by a boundary 
term, they have the same Euler-Lagrange equations. 
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Harmonic maps have been intensively studied, in particular the Dirichlet 
problem for target manifolds with negative curvature has been solved [61] 
|101j |1UU] . and [7T]. The last article is particularly relevant here. However, 
these results do not directly apply in our case because the equations are 
singular at the axis. One of the main technical complications in the proof of 
theorem 3A_ and 3^ is precisely how to handle the singular behavior at the 
axis. 

We present in what follows the main strategy in the proof of the inequal- 
ity (67) (we follow the approach presented in [33] and [35]). The core of the 
proof is the use of a theorem due to Hildebrandt, Kaul and Widman [71] for 
harmonic maps. In that work it is shown that if the domain for the map 
is compact, connected, with nonvoid boundary and the target manifold has 
negative sectional curvature, then minimizers of the harmonic energy with 
Dirichlet boundary conditions exist, are smooth, and satisfy the associated 
Euler-Lagrange equations. That is, harmonic maps are minimizers of the 
harmonic energy for given Dirichlet boundary conditions. Also, solutions of 
the Dirichlet boundary value problem are unique when the target manifold 
has negative sectional curvature. Therefore, we want to use the relation be- 
tween Ai and the harmonic energy Ai' in order to prove that minimizers 
of Ai' are also minimizers of Ai. There are two main difficulties in doing 
this. First, the harmonic energy Ai' is not defined for the functions that 
we are considering if the domain of integration includes the axis. Second, 
we are not dealing with a Dirichlet problem. To overcome these difficulties 
an appropriate compact domain is chosen which do not contain the singu- 
larities. Then a partition function is used to interpolate between extreme 
Kerr initial data outside this domain and general data inside, constructing 
also an auxiliary intermediate region. This solves the two difficulties in the 
sense that now the Dirichlet problem on compact region can be considered, 
and the harmonic energy is well defined for this domain of integration. This 
allows us to show that the mass functional for Kerr data is less than or equal 
to the mass functional for the auxiliary interpolating data. The final step is 
to show that as we increase the compact domain to cover all the auxiliary 
data converges to the mass functional for the original general data. This is 
the subtle and technical part of the proof. 

Finally, we mention that it is possible to construct a heat fiow for equa- 
tions ( 197 )-( 198 ) . Note that the first existence result for harmonic maps used 
a heat flow [53j. In our present setting, an appropriate heat flow that incor- 
porates the singular boundary conditions is constructed as follows. Consider 
functions {a, u) which depend on an extra parameter t. Then, we deflne the 
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following flow 



-2(T 



Aa 



Au-2 



\du\ 



duodrj 
V 



(199) 
(200) 



where a dot denotes partial derivative with respect to t. Equations (199)- 
(200 ) represent the gradient flow of the energy ( 185 ). The important property 
of the flow is that the energy Ai is monotonic under appropriate boundary 
conditions. This flow have been used in |48j as an efficient method for com- 
puting numerically both the solution and the value of the energy at a 
stationary solution. 



5.3 Two surfaces and the mass functional 

In this section we want to define a mass functional over a two-surface S. This 
functional plays a major role in the proofs of the quasi-local inequalities. Let 
us motivate first the definition. 

Consider the mass functional (185) defined over M.^. Assume that on M.^ 



we have a foliation of two-surfaces. For example, take spherical coordinates 



{r,6,(p) and the two-surfaces r = constant. We can split the integral (185) 
into an integral over the two surface and a radial integral. However, the 
integral over the two-surface alone will not have any intrinsic meaning since, 
of course, its integrand depends on r. To obtain an intrinsic expression we 
need somehow to avoid the r dependence taking some kind of limit. For 
extreme Kerr initial data this limit is provided naturally on the cylindrical 
end. The functions a and u have a well defined and non-trivial limit there 
(in contrast with the asymptotic fiat end where they tend to zero). Also, all 
the radial derivatives go to zero at the cylindrical end. This surface define 
an extreme Kerr throat surface (see the discussion in [16] and |15] for more 
details) and it characterized by only one parameter: the angular momentum 
J. Explicitly the functions in that limit are given by 

cro = ln(4|J|)-ln(l + cos2^), coo = -^^^-^ (201) 

1 + cos^ B 

And the area of this two-surface is given by 

Ao = 87r|J|. (202) 
Consider the following functional over a two-sphere 

= j^^ + 4a + smOdO. (203) 
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We can write this integral as an integral over the unit sphere in the 
following way 



1 

2^ 



5-2 



Dar + 4:a + 



\Duj\ 

jj2 



dso, 



(204) 



were dso = sin 6 dddcf) is the area element of the standard metric in and D is 
the covariant derivative with respect to this metric. In complete analogy with 



the discussion in section 5.2, we consider this integral for general functions 
a{9, if) and uj{9, ip) on §^ which are not necessarily axially symmetric. 
The Euler-Lagrange equations of this functional are given by 



Ana 



~9 5 

V 



where 



sin^ 9 e° 



(205) 
(206) 

(207) 



The functional (204) is relevant because the extreme Kerr throat surface 



(201) is a solution of the Euler-Lagrange equations. Equations (205)-(206) 



can be also deduced from the stationary axially symmetric equations (197)- 



(198) taking the limit to the cylindrical end (given by r — )■ in these coordi- 



nates) and using the fall off behavior of the functions at the cylindrical end, 
namely 



a = -2 log(r) + d{9) + 0{r~^), tu = u{9) + 0{r-^). 



(208) 



Note that the 2 that appears in the second term in the left hand side of (205) 
arises from the characteristic — 21og(r) fall off behavior at the cylindrical 
end (we have already discuss this property in equation (189)). Also this 2 



produces the linear term Aa in the mass functional (204). 



The value of the functional (204) at the extreme Kerr throat sphere is 
given by 

M'5 = 8(ln(2|J|) + l). (209) 



The connection between the mass functional (204) and the energy of 
harmonic maps between and is very similar as the one described in the 
previous section for the mass functional A4. Namely, consider the functional 



1 

2^ 



\Dr]\'^ + \Dio\ 



dSr 



(210) 
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defined on some domain Q G S"^, such that Q does not include the poles. 
Integrating by parts and using the identity 



Ao(log(sin^)) = -1, (211) 

where Aq is the Laplacian on we obtain the following relation between 
M and M' 

M'^ = Mfi + A [ log sin edso+ (f (4a + log sin 9) ^^^^^^^^ dl, (212) 
Jn Jan ^■^ 

where s denotes the exterior normal to Q, dl is the line element on the 
boundary dfl and we have used the obvious notation A4q to denote the mass 



functional (204) defined over the domain Q. The difference between Ai'^ and 
A4''^ are the boundary integral plus the second term which is just a numerical 
constant. Note that if we integrate over §^ this constant term is finite 

/ log sin edso = 2 log 2-2. (213) 
Jn 

The boundary terms however diverges at the poles. 

In an analogous way as it was described in the previous section, the 
functional Ai''^ defines an energy for maps [rj, w) : — )■ where denotes 
the hyperbolic plane {(//, w) : fj > 0}, equipped with the negative constant 
curvature metric 

ds^ = dJl±^. (214) 

The Euler-Lagrange equations for the energy A4' are called harmonic maps 
from S"^ — )■ H^. Since A4'^ and A4''^ differ only by a constant and boundary 
terms, they have the same Euler-Lagrange equations. 

The variational problem for the mass functional on the two-surface is very 
similar to the one for the mass functional on described in section 5.2 (see 
[1] for the details). 



6 Open Problems 

In this final section I would like to present the main open problems regarding 
these geometrical inequalities. In the light of the recent results of [72] there 
exists now a very interesting open door to higher dimensions, but this lies 
out of the scope of the present review and hence in this section I will restrict 
myself to four spacetime dimensions. My aim is to present open problems 
which are relevant (and probably involve the discovery of new techniques) 
and at the same time they appear feasible to solve. 
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We begin with the global inequality (54). The two main open problems 
are the following. 



• Remove the maximal condition. 

• Generalization for asymptotic fiat manifolds with multiple ends. 

The situation for the maximal condition in theorems 13 . 1 1 and 13 . 2h esembles 
the strategy of the proof of positive mass theorem by Schoen and Yau |lU2j 
[103] [Tnl] That proof was performed first for maximal initial data and 
then extended for general data. It is conceivable that similar techniques (i.e. 
the use of Jang equation) can be used here also, but it is far from obvious 
how to extend these ideas to the present case. 



The most relevant open problem regarding the global inequality (54) is its 
validity for manifolds with multiple asymptotic ends with non-trivial angular 
momentum. The physical heuristic argument presented in section [2] applies to 
that case and hence there little doubt that the inequality holds. In particular, 



as we already mentioned in section 3.1, for the case of three asymptotic ends 



there are strong numerical evidences for the validity of the inequality 



Theorem 3.3 (proved in [33]) reduces the proof of the inequality to prove 
the bound (57) for the mass functional evaluated at the stationary solution. 
The case of three asymptotic ends (which, roughly speak;ing, is equivalent to 
say that we have two black; holes) is special for the following reason. There 
exists exact stationary axially symmetric solutions of the Einstein equations 
(the so-called double-Kerr-NUT solutions [SI] [22]) which represent two Kerr- 
like black holes. These solutions contain singularities, that prevent them to 
qualify as genuinely equilibrium state for binary black holes. In fact, one 
of the main part in the strategy to prove that the uniqueness theorem hold 
for the binary case is to prove that these solutions are always singular 



However, even if these solutions are singular they can be useful to prove 



the bound (57), because in order to qualify as a stationary point of the 
mass functional Ai all we need is that the functions [rj, u) are regular. It is 
conceivable that some of these exact solutions have this property (that means, 
of course, that other coefficient of the metric are singular) and hence with the 
explicit expression for [t], u) provided by them it will possible to evaluate Ai 



and check the bound (57). In the articles [HZ] [27] the geometrical inequality 



(54) has been studied for these exact solutions. These results provide a guide 



for which subclass of these solutions are potentially useful to prove the bound 



thank Marcus Khuri for pointing this out to me and for relevant discussion on this 
subject 
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(57). Unfortunately the solutions, although explicit, are very complicated 
and it is very difficult to compute the mass functional M. for themj^ 

To compute the value of Ai for this kind of exact solution would be 
certainly a very interesting result which not only will prove the inequality 



(54) for the three asymptotic ends case but also will hopefully provide some 
new interpretation of the double-Kerr-NUT solutions. However this result 
will be confined to the three asymptotic ends case and probably will not yield 
light into the mechanism of the variational problem for the mass functional 
A4 with multiple ends. The basic property which is expected to satisfies 
Ai is the following: if an extra black hole is added, with arbitrary angular 
momentum, then the value of Ai increase. This property is of course another 
way of saying that the force between the black holes is always attractive (in 
particular, it can not be balanced by spin-spin repulsive interaction). The 
variational problem for the mass functional At with multiple ends appears to 
have a remarkably structure. In particular, there is formal similarity between 
this problem and the kind of singular boundary problems for harmonic maps 
studied in [21] . 

We mention in section [2] that there is a clear physical connection between 



the global inequality ( 54 ) and the Penrose inequality with angular momentum 
in axial symmetry. Hence, it appropriate to list here the Penrose inequality 
as a relevant open problem for axially symmetric black holes (for more detail 
on this problem see the review |88]): 

• Prove the Penrose inequality with angular momentum Eq. (17). 



However, it is important to emphasize that it not clear that the techniques 



used to prove theorems 3.1 and 3.2 will help to solve this problem. The 



reason is the following. The proof of the Penrose inequality involve an inner 



boundary, namely the black hole horizon. On the other hand, theorems 3.1 



and 3.2 refer to complete manifolds without inner boundaries. As we mention 



in |3.1| there are results in axial symmetry which includes inner boundaries 



(i.e. [5T] and use similar techniques as in theorems 3.1 and 3.2 (namely, 
the representation of the mass as positive definite integral in axial symmetry, 
see section [5]). However the boundary for the Penrose inequality has a very 
important property: it should be an outer minimal surface (for simplicity, 
we discuss only the Riemannian case). This property is very difficult to 
incorporate in an standard boundary value problem. To see if the mass 
formula in axial symmetry is useful to prove the Penrose inequality with 
angular momentum the natural first step is to prove the Riemannian Penrose 
inequality without angular momentum in axial symmetry using this mass 

thank Piotr Chrusciel for relevant discussion on this point 
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formula. The results presented in [61] [36] contribute in this direction, but 
so far the problem remains open. May be there exists a combination of the 
global flows techniques developed in [7S] [23] for the Riemannian Penrose 
inequality with the mass functional that incorporate the angular momentum 
in axial symmetry. For example, it is suggestive that the strategy for the 



proof of the charged Penrose inequality Eq. ( 15 ) given in [76] and [75] consists 



in first prove the inequality (20) using a flow and a lower bound to the scalar 
curvature that resemble the mass functional. However, a generalization of 
this construction to include the angular momentum is far from obvious. 

We turn now to the quasi-local inequalities. The three main problems are 
the following. 

• Include the charge and the electromagnetic angular momentum in axial 



symmetry for the inequality (64) 



Isoperimetric inequalities in axial symmetry with angular momentum. 



That is, a version for theorem 3.5 (in axial symmetry) with angular 
momentum instead of charge. 



A generalization of the inequality (64) without axial symmetry. 



The inclusion of charge in the inequality (64) is important, of course, 
since charge is the other relevant parameter that characterized the Kerr- 
Newman black hole. But is also relevant for another reason. The angular 
momentum that appears in this inequality is the Komar gravitational angular 
momentum. In the generalization with charge it is expected that the total 
angular momentum (i.e. gravitational plus electromagnetic) appear in this 
inequality. This is shown, including also the magnetic charge, in [56]. In 
addition, this is important in connection with the rigidity statement. There 
is work in progress on this problem [S7] . 



We mention in section [2] that a version of the inequality (64) for isoperi- 
metric surfaces (instead of black hole horizons) could have interesting as- 
trophysical applications, since apparently neutron stars are close to saturate 



this kind of inequalities. Such theorem will be analogous to theorem 3.5 for 



the charge. However, it is by no means clear that similar techniques as the 



one used in the proofs of theorem 3.6 can be applied to that case 



Finally we mention the problem of finding versions of inequality (64) 
without any symmetry assumption. In contrast with the other open prob- 
lems presented here, this is not a well defined mathematical problem since 
there is no unique notion of quasi-local angular momentum in the general 
case. However to explore the scope of the inequality in regions close to axial 
symmetry (in some appropriate sense) can perhaps provide such a notion. 
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From the physical point of view I do not see any reason why this inequahty 
should only hold in axial symmetry. 
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